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ABSTRACT 

Calculations have been made of electromagnetic wave scattering from dielectric disks of 
arbitrary shape and orientation in the high frequency (physical optics) regime. The solution is 
obtained by approximating the fields inside the disk with the fields induced inside an identically 
oriented slab (i.e, infinite parallel planes) with the same thickness and dielectric properties. The 
fields inside the disk excite conduction and polarization currents which are used to calculate the 
scattered fields by integrating the radiation from these sources over the volume of the disk. This 
computation has been executed for observers in the far field of the disk in the case of disks with 
arbitrary orientation and for arbitrary polarization of the incident radiation. The results have been 
expressed in the form of a dyadic scattering amplitude for the disk. The results apply to disks 
whose diameter is large compared to wavelength and whose thickness is small compared to diam- 
eter, but the thickness need not be small compared to wavelength. Examples of the dependence of 
the scattering amplitude on frequency, dielectric properties of the disk and disk orientation are 
presented for disks of circular cross section. 
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HIGH FREQUENCY SCATTERING FROM ARBITRARILY ORIENTED 

DIELECTRIC DISKS 


I. INTRODUCTION 

A solution is presented for tlie scattering from dielectric disks of arbitrary shape and orienta- 
tion. This work was motivated by an investigation of the effects of leaves on microwave remote 
sensing of vegetation (Lang, 1981). Determining the scattering and absorption properties of a 
collection of objects such as leaves requires knowledge of the scattering and absorption properties 
of the individual objects, and in this work the dielectric disk was adopted as a first order model for 
leaves. The disk is also of interest as a model to study the scattering and absorption from other 
particles such as ice crystals in clouds. 

Althougli much work has been done on scattering from perfectly conducting disks (Meixner 
and Andrejewski, 1950; Andrejewski, 1952; Hodge, 1980; Ruck et. al., 1970), very little has been 
done on scattering from dielectric dis*>; and disks with non-circular cross section (Ruck, et. al., 
1970; Bowman et. al., 1969). Howe r, mere has been recent work applicable to thin dielectric 
disks. Weil and Chu ( 1976a-b) developed a numerical appro,Kimation for the resonant region valid 
for thin circular disks. This was done by choosing a set of basisfunctions apropos of a thin circular 
disk and then using a method similar to the moment method to e.xpress the induced currents in 
terms of these basisfunctions. Schiffer and Thielheim (1979) developed a low frequency approxi- 
mation for the thin dielectric disk by using a Rayleigh type approximation. In this solution the 
induced currents are found by assuming that the internal fields have the same form as would be 
obtained in statics. The approximation is valid for disks whose physical cross section is much larger 
than the thickness and whose thickness is much less than a wavelength. The solution applies at 
high frequencies (i.e. to disks whose cross section is large compared to wavelength) only if the disk 
is sufficiently thin. 

The work to be described here applies specifically to disks whose cross section is large com- 
pared to wavelength but is not restricted to disks which are thin compared to wavelength nor to 
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disks of circular cross section. The solution is obtained by using a variation of the Kirchoff approxi- 
mation employed in physical optics scattering from rough surfaces. In this approximation the fields 
on the surfaogi are approximated by the fields on a plane tangent to the surface. In the application 
to the disk, the fields inside the disk have been approximated by the fields inside a slab of the same 
orientation and thickness. Given the fields inside the disk the scattered fields can be obtained from 
currents induced by these internal fields. In the required integration the cross section of the disk 
need not be circular nor must the disk be thin compared to wavelength. 

In the sections to follow a formal solution for the scattered electric fields is developed in terms 
of the unknown fields inside the disk. Then the Kirchoff-style approximation for the internal fields 
is derived for an arbitrarily oriented disk. Using this approximation the volume integration is 
carried out to obtain explicit forms for the scattered fields for observers in the far field of the disk. 
Finally, this solution is expressed in the form of a dyadic scattering amplitude for the disk. Ex- 
amples of the scattering amplitude are presented to illustrate the effects of frequency, dielectric 
properties and orientation of the disk. Several checks have been made on the solution and are also 
described. It can be shown that the solution satisfies conservation of energy at high frequencies and 
that the solution predicts the accepted result for radar cross sections in the special case of normally 
incident plane waves. 


II. GEOMETRY 


The problem to be addressed here is to calculate the fields scattered by a plane wave of 
arbitrary polarization incident on an arbitrarily oriented disk. The problem is illustrated in Figure 
1. The plane wave is assumed to have polarization (direction of the electric field) q and to be prop- 
agating in the i direction. Specifying i in terms of the spherical coordinates and for con- 
venience letting the incident wave have unit amplitude, one can write: 

Ejnc q exp(j ko 1 • r) (I) 

where 

I = - sin(0j) [cos (0j) X + sin(0j) y] - cos(0j) z (2a) 

and kg = Aside from the requirement q • i = 0 the unit vector q is arbitrary. The orientation 

of the disk is described by the Eulerian angles (d,<f>,y) as indicated in Figure 2. (See Goldstein, 1 966 
and also Appendix B). Starting with the axes of the disk (x', y', z') aligned with the reference 
system (x,y,z) the disk is rotated to its arbitrary position by: 1) rotating ^ degrees about the z' * 
z axis; 2) rotating about the new x'-axis d degrees; and 3) rotating about the new z' axis y degrees. 
All rotations are counter clockwise. In terms of these angles the normal to the disk, n, is: 

n * sin 6(sin0 x - cos 0y] + cos 6 z (2b) 

The shape of the disk is described by the function S(r'), defined in the coordinate system (x', y' z') 
fixed on the disk. It is the cross sectional shape looking along the local normal to the disk (i.e. 
z' axis). The disk is imagined to be cut from two parallel planes T meters apart by a cookie cutter 
with shape S(x',y'). S(x',y') = 1 on the disk and zero otherwise. 

The problem as stated above is sufficiently general to describe scattering from a plane wave of 
arbitrary polarization and direction of propagation incident on an arbitrarily oriented disk. The 
orientation of the disk, the shape of the disk, S(x',y'), its thickness, T, and its dielectric properties 
are all arbitrary in the analysis to follow. 
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III. THE FORMAL SOLUTION 

To obtain a solution for the scattered fields in terms of the fields inside the disk the constitu* 
tive relations D = e^E + P and B * Mq + M) are first used to obtain a form of Maxwell’s equations 
with P, M and conduction current, as sources. Then, Fourier transforming and introducing the 
magnetic vector potential, A, one obtains the wave equation: 


A + k^A * - - jcoF + V X M) 


(3) 


which has the particular solution 


ACT)=^ j’tJ„-juP + rXM)^^dr' (4) 

where R « 1 7 - r' | and dF * dx'dy'dz' and the integration is over the volume of the disk. The 
scattered electric field can be obtained from the vector potential in the form: 


E(f') = jw(A + ~^(V- A)] 
k2 

which yields (Substituting Equation 4 into Equation 5): 


(5) 


E(r')=jw|ao j‘jT(r')-G(?/r')df' 

V 

In Equation 6, TjCf) =T^ -jwF + V X M and S(T/T') is the dyadic Green’s function for free space: 

k 


( 6 ) 






_ I ^kgR 


(7) 


4ffR 


To complete the solution the sources jw P and V X M are required. In principle, arbitrary 
relationships between Jg, P and M and the fields E and B can be used; However, for most media 
important in remote sensing of the earth (e.g. ice, snow, water, vegetation, soil, etc.) it is appro- 
priate to assume a simple linear relationship: 
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( 8 ) 


T, - aE 

P = e„XeE 

M - XmB 

where is the electric susceptibility and x^j the magnetic susceptibility. In most cases Xm ^ 0. 
Assuming this to be the case the total current Tj becomes 

Xj. = -jw6o E (9) 

where Xg is an equivalent, complex, electric suspectibility: 

Xe*Xe+jo/wCo (10) 

The complex relative permittivity of the medium, is: 

r, = Xe+1 01) 

* €j+ja/oj€^ ( 12 ) 

where is the relative permittivity of the disk when a = 0. Substituting Equation 9 into Equation 
6, one obtains an expression for the scattered electric field in terms of the electric field E(t 0 inside 
the slab: 

^scat = 1^5 C ^ (7') * ^ O'/T') dT' (13) 
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IV. THE FIELD INSIDE THE DISK 

Equation 13 applies for observation points both inside and outside of the object. If the 
observer is inside the object, then Equation 13 is an integral equation for¥(r). When the observer 
is outside, then Equation 13 is a solution for the fields scattered by the object in terms of the 
fields inside the object. The problem, of course, is that the Helds inside the object are not known in 
either case, and appro.ximations are necessary. The approximation to be used here will be to 
assume that the Helds inside the disk are the same as would exist in a slab of inifinite parallel faces 
of the same thickness and orientation. This is an approximation which ought to be reasonable for 
disks whose cross section is large compared to wavelength and thickness and is independent of 
the dielectric constant of the medium. It is the analogue of the tangent plane (Kirchoff) approxi- 
mation employed in the physical optics approach to scattering from surfaces. 

The fields inside the equivalent slab are most easily obtained in the coordinate system (x V'z') 
fixed on the disk. Once obtained they can then be expressed in terms of the reference coordinates 
(x,y,z) by a suitable transformation of coordinates. The problem as seen in the coordinate system 
fixed on the disk is shown in Figure 3. Notice that the origin of the primed coordinate system is 
at the center of the slab and that the plane wave is incident from the arbitrary direction, i . As 
seen by an observer on the slab the incident plane wave has the form: 

where primes (') are being used to remind the reader that these are coordinates in the reference 
frame of the disk. To find the fields inside the slab, it is convenient to resolve the incident wave 
into its horizontally and venicaJly polarized components as seen by an observer on the slab. Letting 
h and v be the appropriate unit vectors in the slab’s reference system, one can write q in terms of 
its projections on these polarization vectors: q = (q 'hjli +(q *v)v. With this notation a unit 
amplitude incident wave, as seen by the observer on the slab, can be written: 

^inc *h)h + (q • v )v ] (15) 
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The fields inside the slab are also plane waves. Explicitly separating their dependence on the co< 
ordinate perpendicular to the slab (z') and coordinate parallel to the slab boundaries (x',y'), they 
can be written in the form: 

^slab (f') * (x'y)e-J(nkoK;)z' ( 16 a) 

where 

fi^Xx^yO “ e+ ^ Ky y I (16b) 


n is the index of refraction of the medium in the slab (n -y/e^) and K = x' + Ky y' + z' 
is a unit vector in the direction of propagation of the wave inside the slab wliich propagates in the 
+ z-direbtion (i.e. > 0). The vector amplitudes e+ of the waves inside the slab can also be re- 

solved into horizontally and vertically polarized components. Thus> let h| , v|, be unit vectors 
in the direction of the electric field inside the slab when the slab is excited by an incident hori- 
zontally or vertically polarized wave, respectively. (For example, if q = h then the waves inside 
the slab will have polarization h* .) Also let the (scalar) amplitude of the waves inside the slab 
due to an incident horizontally or vertically polarized wave of unit amplitude be e^, e^ respec- 
tively. Then, in the case of an incident wave of arbitrary polarization q one can write the ampli- 
tude of the waves inside the slab as , 

e± * (q • h ) ej ii* + (q • V ) e; V* (17) 

The scalars e^ ^ are characteristics of the slab, independent of its orientation. They are obtained 
from the boundary conditions at the slab interfaces in the same manner as is conventional in find- 
ing the reflection or transmission coefficients (e.g. Bom and Wolf, 1959). One obtains: 


+ ^ h,v ^h,v ^ 


(18a) 



‘h.v 






(18b) 
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where 




2ko T n k; 


(1 9a) 


A* 


~k,TK±nK;j 


(19b) 


I “ Xl, V 


(1 9c) 


'i+Vv y\ 



(19d) 


-n^ 


(I9e) 


1 K 

\ 

^ n ii 


(190 


where = K • h > 0 and = I * h < 0; and the rj^ ^ and tjj y are the reflection and transmission 
coefficients, respectively, of a half-space with the same dielectric properties as the disk. The minus 
sign (-> in Equations 19e, f are the result of the convention adopted here for i and h (i.e. the wave 
is incident on the slab from the h side). 
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V. THE SCATTERED FIELD 

The fields scattered from the disk can now be obtained by substituting the fields inside the 
disk (Equations 16) into Equation 13. To do the integration in Equation 13 over the volume of 
the disk it will be assumed that the observer is many wavelengths from the disk (k^R » 1) and 
that the disk is small compared to the distance between the disk and the observer (e.g. L/R « 1 
and kL~/R « 2v where L is a dimension characteristic of the size of the disk). In this case one 
can approximate the Green’s function by; 




( 20 ) 


where R^ is the distance from the center of the disk to the observer and 6 is a unit vector from the 
origin to the observer (i.e. in thl f.'.irection of propagation of the wave scattered to the observer). 
Substituting Equation 16 into Equation 13, using Equation 20, and separately doing the integra- 
tion over the (x', y') and z' coordinates one obtains the following compression for the fields scatter- 
ed from homogeneous disks (Xe and a constant): 


— T — 

^scat “ ^0 


R(T)S(P;) 


( 21 ) 


where 


T/2 

R(T)* J 6 x[[g+6j*'o[nKi-oi]z' + ^_5}ko[nK;+oilz'jx6 jdz’ 
-T/2 


= T { 6 X [5+ X 6 J sine ['/sk^T (nK^ - op] + 
+ 6 X [e, X 0 ] sine [V7k^ T (nKj + op]} 

oe 

S(^;) * y y S(rpej2ff 


(22a) 

(22b) 

(22c) 


The subscript “t” in Equation 22c denotes transverse components; r[ = x' x' + y' y' 
k — k 

and v[ * [nK^ -op*~h X[nK-6)Xh]. Notice that using Snell’s law at the slab 
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(22d) 


boundaries, n K| f!, so that one can write V[ in terms of i and o : 

((I-6)Xh| 

^ 2ir 

The vector direction of the scattered field, is determined by the vectors 6 X (e+ X 6 ) 
associated with 1^(T). Since 6 is a unit vector in the direction of propagation of the scattered field 
(i.e, pointing from the disk to the observer), the scattered fields are clearly transverse as seen by 
the observer. This is as expected since a far field limit has been taken to obtain the solution. 


It is conventional to express the scattered Helds in terms of their vertically and horizontally 
polarized components. To do so, let p be a unit vector in the direction of vertical or horizontal 
polarization (defined by the observer). Then p scattered electric field of polari- 

zation p due to an incident wave of polarization q (as defined at the transmitting antenna). Noting 
that p • (6 X (5^ X 6 )] * p • e+ and using Equation 17 for e+ one can write: 

ik R 

P-E,„,(4) = (P (23) 

where the dyadic F(i ,6 ) has the form: 


F(1 »o ) = (hgh + v7 V e”I sine [Yzk^ T (nK' - opj + 

+ [hgh ej + vjv ej] sine [Yzk^TinK!^ - op] 
Equation 23 suggests the definition of a dyadic scattering amplitude f(i ,6 ): 


(24) 


f(i‘,6) = (25) 

47T 

To complete the solution it is necessary to find explicit expressions for the polarization 
vectors v*, h| and v, li associated with the slab and with q, p associated with the incident and 
scattered waves. To do so the following definitions will be adopted: A horizontally polarized 
wave (h) is one whose electric field is perpendicular to the plane of incidence (plane defined by the 
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z-axis and the direction of propagation) and a vertically polarized wave (v) is one whose electric 
field is parallel to the plane of incidence. In addition, v, li, k form a right hand orthogonal set 
where k is a unit vector in the direction of propagation of the waves: v X h - k. With these 
definitions, one may write explicit forms for the polarization vectors in the reference (Toofdinate 
system. Thus, for the incident wave (i.e. q * hj or Vj) one has: 

r ^ I X z 

"i ,777. (■:. 


Vj * lij X I (26b) 

where i is defined in terms of the spherical coordinates (6^, 0j): * 

i * - sin 6j(cos x + sin y] ~ cos 6j z (26c) 

and in the case of the scattered fields the unit vectors hj, Vj in the direction of horizontal and 
vertical polarization (q * hj or v^) are: 


- 6 X z 

h ~ ■ 

(27a) 

* 1 6 X z I 


Vj^^hjXd 

(27b) 

where the spherical coordinates (Sj.^j) have been adopted to define 6: 


6 = sin 6Jcos x + sin y] + cos 6^ z 

(27c) 


The unit vectors h, v are the direction of horizontal and vertical polarization as defined by an 
observer on the slab but expressed in terms of the reference coordinate system (unprimed). Since 
h = z' as seen in the reference system and since the direction of propagation of the incident wave 
as seen in the reference frame is 1, one obtains: 


1 X h 
liXhl 


(28a) 


V * h X i 


(28b) 


where the normal to the slab h is defined by the Eulerian angles (6,(p,y): 

n*sin6 [sin <^x-cos0y] +cos^ z (28c) 
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The unit vectors v* for the direction of electric field inside the slab are also required. Letting 
K+ be a unit vector in the direction of propagation of the two waves inside the slab one has: 



j * K* X h 
*'«*lKiXh| 

(29a) 

1 


(29b) 

It is shown in Appendix A that 

*• 

+ a] 

n 

(30a) 

where 


s= - (i • h) ±/ivr-\) + (i • n)‘ 

(30b) 

Thus one obtains: 

. 



hf*h 

(31a) 


v*»“{v + n+ h X n] *- [v + J2x] 

C n * n 

(31b) 


With the polarization unit vectors given above one can write the solution for the scattered , 

field in the final form: 

« 

P • = tP • “(i, 6) • ^32a) 

T(1, 0) « S (»;) P (1 6) (32b) 

4JT 

F(i, o) = [h h + - (v + n«) v e^] sine [‘/zko T(nK 2 + o^)! 

" 32c) 

+ [h h ej + i (V + n+) v e +] sine [Vzk^ T(nK; - o^)] 
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I 


where the e^ y are given by Equations 18 and 19 and 

i' = I • h * - sin sin 6 sin (<p - 0j) - cos cos 6 

0*2 * 6 • h * sin 6^ sin 0 sin (0 - ^j) + cos 0^ cos 0 
nK^ * nK • h =y/(n^ - 1) + (i • h)“ 

= Hi ± nKi) h X n 

Explicit forms for the shape function S(Pj) are £^ven in Appendix C. 


(33a) 

(33b) 

(33c) 

(33d) 


ii 
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VI. EXAMPLES 

Examples of the scattering amplitude will be presented in this section for a sample disk to 
illustrate some of the characteristics of Equations 32. Calculations have been made to illustrate 
effects of frequency, disk orientation and dielectric constant on the scattered fields. The calcula- 
tions are for a nominal disk of radius 10 cm and thickness 0.5 cm, and except where indicated the 
frequency of the incident radiation is 9 GHz and the dielectric constant of the disk is = 25 + j 1 1 
which is representative of leaves in the microwave region. In these examples the incident radiation 
has been choosen so that the direction of incidence is In the y-z plane 30° from the z-axis (i = 
■“cos(30°)z sin(30°)y) and the normal to the disk, n, is in the plane of incidence (n = -sin0 y + 
cos4> t). The geometry is shown in Figure 4. To do the calculations the observer is imagined to be 
on the surface of a large sphere centered on the origin and the scattering amplitude p*f(t,6)*q 
is computed as the- observer is moved around one of the great circles formed by passing a 
plane through the z-axis. The magnitude of the scattering amplitude, | p*f(I,6)*q |, has been dis- 
played as a polar plot of amplitude versus angle on the great circle (Figures 4- 15). In each case the top 
of the figure corresponds to the point when the observer is on the z-axis. When the observer is in 
the plane of incidence (e.g. Figures 6-15) backscatter occurs at 30° to the right of the top and 
forward scatter at 150° to the left. The polar plots are all linear in amplitude. The outer circle 
represents | p*^i,6) *q | = 1.0 and the inner circle is | p*T(i,6) *q | = 0.2. This scale is the same in all 
the examples to be presented here (Figures 5-15) except for the 45° and 90° cuts in Figure 5 where 
the outer circle vas chosen to be | p*?^i,6)’ q i = .05 to better show the detail of these small 
scattered fields. 

The elements of the scattering amplitude most commonly of interest are: 

k • * hj = fhh 

Vj • 7(f,o) • hj = fyh 

hj • 7(1,6) • Vj * fj^y 
Vj ‘ 7(1,6) • Vj = fyy 


14 


The elements ^hh and are proportional to the horizontal and vertical polarized scattered field, 
respectively (as seen by the observer), due to an incident wave of horizontal polarization. The 
elements, ^hv and fyy are proportional to the scattered field oi‘ horizontal and vertical polarization, 
respectively, due to a vertically polarized incident wave. The unit vectors, hj, Vj and hj, Vj are 
indicated in Figure 4 for reference. The magnitude of these four components of the scattering 
amplitude have been calculated. 

Figure 5 shows the four components of the scattering amplitude fj^j^, fy|^, fj^^ and f^^ for an 
observer in the plane of incidence (0“), in a plane perpendicular to the plane of incidence (90®) 
and in a plane half way between these two (45®). (See Figure 4.) The calculations are for the case 
n = z. Each plot indicates the magnitude of p*f(f,o)*q in polar form with the z-axis at the top 
of the polar plot. Notice that ^vh * ^hv ® when the observer and n are in the plane of incidence 
(0®). That is, there is no depolarization in this case. This is true whenever 6, h and i are coplanar. 
Also notice that ^hh “ ° when the observer is in a plane perpendicular to the plane of incidence. 

This is a coincidence dependent on the definitions choosen for the polarization unit vectors. In the 
plane of incidence (0®), the scattering amplitude has two major peaks, one in the direction of inci- 
dence, 1 , (the downward pointing peaks in Figure 5) and one in the specular direction (the upward 
pointing peaks). The specular peak is in the direction of radiation reflected from an identically 
oriented infinite slab and the downward or “forward scattered’’ peak is in the direction of radiation 
transmitted through such a slab. For a highly conducting disk, the forward scattered radiation com^^ 
bines with the incident field to produce a “shadow” behind the disk. 

In the physical optics regime the shape and amplitude of the forward scatter and specular 
peaks depend on the product, ka = 27ra/X where a is the radius of the disk and X is the wavelength 
of the incident radiation in the ambient medium (Appendix C), This is illustrated in Figures 6-7 
which show ^hh and fy^ (respectively) for several values of ka. The geometry is as illustrated in 
Figure 4 with the observer in the plane of incidence (0® plane). Notice that as ka increase the peaks 
become larger and narrower in the directions of the reflected and transmitted waves. For a circular 
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disk, the peaks behave like Jj(x)/x where J^Cx) is a Bessel function of first kind and x depends on 
ka and disk orientation (Appendix C). 

Changing the orientation of the disk affects the direction and shape of the peaks. This is 
illustrated in Figures 8 (fj^|^) and 9 (fyy) for the special case where o,n i, z are all coplanar. In 
these examples the observer is in the 0^ plane (Figure 4) and h is also in the O’’ plane but at several 
different angles with respect to the incident wave. The four polar plots have been obtained with h 
parallel to i (5 = 0*) and then moved to the left so that the angle between h and i is 15*, 45* and 
finally 60*. Notice that as the direction of h is changed the specular peak changes direction but 
that a peak always remains in the forward scatter direction. As the disk appears more edge on to 
the incident radiation (9 > 60*) the specular and forward scatter peaks merge and their amplitude 
decreases toward zero. (The theory is not applicable near grazing incidence.) 

The effect of the dielectric constant of the disk on the scattered fields is illustrated in Figures 
10-15. The geometry is as illustrated in Figure 4 with the observer in the 0* plane and h = z. The 
frequency is 9 GHz (ka = 2). Figures 10 and 1 1 are for a lossless dielectric (?J, real) and show the 
effect on ^hh and fyy, respectively, of increasing the relative dielectric constant of the disk. With the 
choosen frequency (9 GHz) and thickness (0.5 cm) the disk goes througli a quarter wave resonance 
near ?^ = 11. These figures show the effects on the scattering amplitude of increasing from 1.0 
(upper left) through resonance and then to large, non-resonant values. Initially the scattered 
radiation increases as increases; however, near the resonant condition, e j. = 1 1, the reflected 
radiation decreases markedly. The forward scatter peak, on the other hand, is only slightly affected 
and appears to increase somewhat near resonance. This pattern of decreasing specular peak occurs 
at all quarter wave resonances, which for this disk (T = 0.5 cm) and frequency (9 GHz) occur at Cj. = 
1 1, 44, 99, etc. For large values of Ty not near resonance the scattering amplitude approaches a 
limiting shape as illustrated by the last two examples in Figures 10 and 1 1 (?j, = 70 and 1000). 
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The presence of loss in the dielectric significantly affects the behaviour of the scattered fields 
near resonance. With even small amounts of loss, the resonance is virtually eliminated. This is 
illustrated in Figures 12 and 13 which show examples of the scattering amplitude ^hh and fyy, 
respectively for the nominal disk with fixed loss (Im - 2) and with the real part of having 
values near the quarter wave resonance at e'y * 1 1. For values of Im less then 2 a decrease in the 
specular peak is noticeable and for values of Im greater then 2 the resonance is even less 
apparant than in Figures 12 and 13. 

Figures 14 and 15 present examples of the scattering amplitude for the nominal disk with a 
purely lossy dielectric: Re ~ 0. Four examples illustrating the effects of increasing Im 6"^ = e" are 
shown. The behaviour is very much like increasing Re except for the lack of resonances. The 
scattering amplitude quickly reaches a limiting value and is essentially independent of for 
10 . 
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VII. NORMAL INCIDENCE 

The special case in which the radiation is normally incident on the disk provides a check on 
the solution. In this case one expects the horizontally polarized radiation scattered from horizon* 
tally polarized incident radiation to equal the vertically polarized radiation scattered from vertically 
polarized incident radiation, and that there will be no cross polarized scattered radiation. That is, 
one expects: 

^hh “ ^vv 
fhv 

Also, at normal incidence one expects the backscatter cross section, to be 

®N * R; 

where is the backscatter cross section of a perfectly conducting disk and is the reflection 
coefficient of a slab of the same thickness as the disk (Ruck et al., 1 970). It will be shown that 
these results are obtained from the solutions developed here. 

At normal incidence n * z and h X n = v, and consequently, the dyadic F(I,o) has the simple 
form: 

. ^(i,o) * [ej^ sincfd'*’) + ejj sinc(d“)J + ^ [e^ sincf^'*) - e“ sinc(0‘*’)] (34a) 

where 

« >/2ko T [nKj ± o; j * «/2ko T ( n ± 1 ] (34b) 

The polarization vectors at normal incidence are hj » hj = h and - Vj * + Vj » + v where the minus 
sign (—) in the vertical polarized scattered radiation v^ is a consequence of the convention chosen 
for polarization (Vj * hj X 6) and the fact that a normal incidence *i * - 6. With these polarization 
vectors one has: 

itj • F • hj * ej^ sine id*) + ej^ sinc(5“) 

hj • F • Vj » 0 

(35) 

Vj • F • hj = 0 

Vj • F • V- = e“ sinefd''’) - sinc(6”) 
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Clearly * fyj^ * 0 at normal incidence. To determuiw* thw like>polarized components of scat- 
tered electric field, ^hh and fyy, note that at normal incidence 


'h,v— 

(36a) 

1 - n 

(36b) 

and also that the reflection coefficients R|^ y of a slab of thickness T which in the notation adopted 

here are*. 

- 

^h,v 

(37) 

have the following form at normal incidence 


Rh - y 2 lU + n) ej ej“" + (1 --n) ej; e"j« 

(38a) 

Ry = «/ 2 [(l + n) ej ei“" - (1 - n) e^ e"J“"] e“i® 

(38b) 

where a = '/zk^T. Now substituting Equations 36 into Equations 35 and using Equations 38 one 

obtains: 



(39a) 


(39b) 


And since Rj^ * - Ry * Rjj it is clear that * - fyy. The sign difference between f|^^ snd fyy 
is a consequence of the convention chosen for polarization: vertical polarization for scattered and 
incident waves along the same line of sight are in opposite directions. 


The radar cross section at normal incidence is readily obtained from the preceeding results 
and also provides a check on the solution. Recalling that the definition of radar cross section at 
normal incidence is 
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one obtains 


( 41 ) 

Now using the result that n^ - 1 and that at normal incidence S (F^) = where is the 
area of the disk (Appendix C) one obtains: 


if2 c2 
IT 

The radar cross section of a perfectly conducting disk at normal incidence is 
the radar cross section of the dielectric disk at normal incidence can be written 

®N*IRn 

which agrees with the expected result (Ruck, et al., 1971). 


(42) 

k^ Sj/ff and so 

(43) 


This result is especially simple for thin 

disks. (k(jTn « 1). Keeping only lowest terms in k^T in Equation 38 one can show that 

« 

RNaj!^(n2-l) (44) 

ar.U so at normal incidence 

( 45 ) 

That is, the thin disk ai normal incidence is Rayleigh-like, scattering power proportional to 
rather that than as 1/X^ as the perfectly conducting disk would suggest. 



VIII. ENERGY CHECK 


A check has been made to see if the solution for the scattered fields conserves energy. This 
has been done by showing that the solution satisfies the “optical theorem”. The optical theorem 
relates the total power scattered (Pj) and absorbed (Pj^) to the scattering amplitude in the forward 
direction (Bom and Wolf, 1959). One can state the optical theorem as follows: 

where Pj is the power density of the incident wave, i is a unit vector in the direction of propagation 
of the incident .vave, and q is a unit vector in the direction of polarization of the electric field of 
the incident wave. Defining scattering and absorption cross sections by = Pj/Pj and - P^/Pj 
respectively, the optical theorem becomes: 


where 


<7, • q) • [f(i,6) • q] ’ dS2 

o..i- 


' 1 Sphere 


P; 



a E • E* dv 


(48a) 

(48b) 


disk 


the integration in the expression for is over solid angle df2 = sind ddd^ and the fields E in the 
expression for are the fields inside the disk produced by a unit amplitude incident plane wave. 
They are given in the reference frame fixed on the disk by Equations 16 in the text. For purposes 
of doing the energy check, the orientation of the disk is unimportant and can be chosen for con- 
venience; the choice in which the primed and unprimed systems are identical is most convenient 
(6 = 0 * 7 = o). In this case, using Equations 16 in Equation 48b one obtains: 


j[l ^ I' + I'] 2 R, [e- ej* (qj.q;*) j (49) 


where Sg is the cross sectional area of the disk and q e jh,v} . The are given in Equations 18a,b, 
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and Ky * Re(nkgKp, * ImCnk^ K^) where is given by Equation 33c, 

An analytical and a numerical method has been used to check that the solution given by 
Equation 25 does satisfy the optical theorem when ka»l, The analytical check is made by 
asymptotically evaluating the integral for a, given by Equation 48a for ka»l . This evaluation is 
carried out in Appendix D where an explicit formula for is given by Equation 6D. When this 
expression is used in the optical theorem along with Equation 49 and Equation 25, the conservation 
theorem has been shown to be satisfied identically. 

Following the analytical check, a numerical check was made by computing and from 
Equations 48 and 49 and comparing their sum with the right hand side of Equation 47. Calculation 
were made for the case of a circular disk of radius 7cm, thickness 1mm and relative dielectric 
constant ej. * 36 + jl3. The dielectric constant of the disk was choosen using the de Loor formula 
(Fung and Ulaby, 1 978; de Loor, 1968) for leaves with 70% water at a frequency of 7 GHz. Two 
expressions for the total cross section were developed: The first, Oj, was found by computing 
and Oj from Equation 48 and 49. The second, Cj, was obtained by taking the imaginary part of 
4 jt (q*f(i,i)*q)/kQ. If the optical theorem is satisfied exactly then aj « Oy. The results of these 
calculations are shown in Tables I and II. Table I shows the results for a horizontally polarized 
wave incident at 30** with respect to the normal to the disk and Table II shows the results for a 
vertically polarized wave also incident at 30®. The tables list a^,a^,Oj = and the albedo, 

o^/oj all obtained numerically from Equations 48 and Oj obtained from Equations 47. 

Examples are shown in each table for frequencies of 1 GHz, 4 GHz and 7 GHz corresponding to ka 
of about 1.5, 5.0 and 10.3 respectively, since the theory is a high frequency theory valid for large 
ka, one expects energy conservation to be satisfied best for large ka. The results indicate agreement 
which improves with increasing ka and is within about 2% at the largest value of ka. To obtain 
the numerical results for scattering cross section the integrand, | f (l,o) * q |" in the expression for 
Cj has to be computed on a unir sphere. The results shown in Tables I and II were the result of 
computations using a grid of equally spaced points separated about 3® in d and 0. 
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APPENDIX A 

Polarization Vectors inside the Slab 

The propagation vector for waves inside the slab can He obtained from the boundary 
conditions (SnelTs law) at the surface of the disk. To do so, adopt the following notation in the 
reference frame of the disk: 

Rt - n ko [K;f x' + z'J * n ko K * A1 

7«ko li^x' + iy y' + i^i'J *ko 1 A2 

Since the boundary conditions must be satisfied for all x', y \ one has: 

i; • n k„ Kf i n k„ Ki A3 

k„i;,-nk„K*'4 nkoK^, A4 

and since by definition the waves in the slab are propagating in opposite directions along the z'-axis, 
it follows that; 

R±*nko [K;x' + K;,y'±K;z'] A5 

Now using K+*K+ * (n kg)“ which follows form the wave equation, one can solve for K^: 

y(nk„)2 [I-(K;)2-(K')2] 

n Kq 

- v^!lv^-kif(i;)2+(ip2j 

n Kq 

=■ -v/(n^-l) + (>i)^ 

n * 

* -y/(n^ “• 1) + (1 * h )• A6 

n 

Now using Equations A3, A4 and A6 in Equation AS one can obtain an expression for’K^: 

* kg [ix x' + iy y!J ± n k^ (n- - 1) + (i * h jn 

= k^ I n - (i * h ) n I ±\/(n-- I) + (i • h)^ h } 

«k^{ 1 h } A7 

where 

n+*-(i • h ) ±>v/(n- ~ I) + (i • h AS 
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Thus 


Kj. *- [l + fli nl A9 

This result can now be used to obtain an expression for the unit vectors in the direction of hori- 
zontal and vertical polarization inside the disk. One obtains: 


k+Xn iXft , . 

ft* X ■ X - . ^ ft 

* |iq.Xn| liXnl = 

V* = hi X Kx * h X kx = i [ V + h X h J 

e e X I 


AlO 

All 


Notice that although h^ * h”= h the unit vectors of the vertically polarized waves in the slab in 
general are not in the same direction: v^ # v”. This is true because the two waves inside the slab do 
i.ot propagate in the same direction. This is illust.'ated in Figure 16. At normal incidence = - 
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APPENDIX B 

Euler Angles for Orientation of the Disk* 


The Euler angles which specify the orientation of the disk consist of three coordinate rota* 
tions. Starting with the reference coordinate system and coordinate system for the disk aligned, 
one performs the follosving counterclockwise rotations: 
t . Rotate ^ degrees about the original z» a' axis 


“x’ “ 


" COS0 

sin0 

0 “ 


"x 

y’ 


-^in0 

COS0 

0 


y 

-z' - 


0 

0 

1 


-2 . 


2. Next rotate d degrees about the new X-axis 


“x"“ 


“ i 

0 

0 " 


“x'“ 

y" 

ss 

0 

cosd 

sinO 


y' 



0 

sind 

cos^L 


-Z'- 


3 . Finally rotate 7 degrees about the new 2-axis 


- Itf- 


co$7 

sin7 

0 


”x"" 

y»" 

m 

~sin 7 

COS7 

0 


y'' 

.z'"- 


- 0 

0 

1 « 


-z". 


These three rotations are applied in succession to obtain the relationship between the reference co- 
ordinate system (x,y,z) and the final coordinate system fixed on the disk (x'.y'.z'). Thus: 


” x' " 


COS7 

sin7 

0“ 


“ 1 

0 

0 ” 


0 

0 

1 

siji0 

0 " 

y' 

3E 

--sin 7 

COS7 

0 


0 

COS0 

sin0 


*“Sin0 

COS0 

0 



- 0 

0 

1 - 


.0 

“Sin6 

COS0 . 


. 0 

0 

1 „ 


x' * tcos7cos<^-sin7sin0cos^l (cos7sin0+sin7cos0cos01y + {sin<9sin7]2 
y' - Isin7cos0+cos7sin0cos^l x + (““SinTsin^+cosTCos^cos^^ly + [sin 0 cos 7)2 
2' * (sini 9 sin 0 ) x + (--sin0cos0] y + ' . )$0] 2 
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APPENDIX C 


The Shape Function S(v[) 

The shape function is defined in the reference frame of the slab by 

§(:p;) = J J s(f;)e"i2’^^l -r; 

where S{rp is the cross sectional shape of the disk and where 

7| * x' X + y' y 

It is necessary to express this function in the coordinates of the reference frame (but the integration 
is done in the primed coordinate system.) The problem is to express in terms of the reference 
coordinates (angles 6^, and d^, and 6, <f>, y). Since and Vy are the projections of on the 
x' and y' axes and since dot (scalar) products are invariant under coordinate rotations, one can write 

2ir 

where i, 6 are the directions of propagation of the incident and scattered waves as seen in the 
reference coordinate frames; 

i = ~ sin0j [cos0j x + sin^| y] - cos z 
6 = sindj [cos <^5 X + sin^j y] + cos0j z 

and x', y' are unit vectors along the x and y axes of the disk but as seen by an observer in the 
reference system. These unit vectors can be obtained from the Euler angle rotations (Appendix B): 

x' * (cos7 COS0 “ sin y sin <p cos <?] x + 

+ (cos 7 sin0 + sin y cos0 cos5] y + sind sin 7 z 
y' = “ Isin 7 cos0 + cos 7 sin<f> cosd] x + 

+ sin 7 sin<^ + cos 7 cos<^ cos5] y + sin5 cos 7 z 
z' = sin^ [sin0 x - cos^ y] + cosd 2 
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from which one obtains 


where 


*'x 2ir ^ 

' X 

*'y * 2ir IP sinT + « cosy] 

a * cos(?[sindj sin (^ 0 j) + sin 0 j sin (^^ 5 )! -sind [cos^l + cosdj] 
0 = sin (9j cos + sin^j cos ( 0 -^j) 


An important special case is that of a circular disk of radius a. By changing variables $0 that 
the integration can be done in cylindrical coordinates one obtains: 


where 


S'O'I) = (2Jr»'a) 


v = ^ vA" + 
2ir^ 


and J j (x) is a Bessel function of first kind of order unity. 
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APPENDIX D 


High Frequency Evaluation of cTj 


The scattering cross section, is evaluated asymptotically for large k^a in this appendix. 
Starting with the definition of scattering cross section (Equation 48a) and using Equation 25 for 
f(i,o) on obtains: 


a 


s 


C y |S(F;)|2 w(‘»)(0,0)dJ2 
sphere 


1-D 


where 


W(4) (a, 0 ) = I F(i,o)*q P = I sine + sine (0")c"q" i" 2-D 


0* = ‘/2kJlnK;±o;i 


3-D 


c 


= 




4-D 


and q e | h.v} and the subscript “e” in Equation 2-D indicates polarization unit vectors of waves 
inside the disk (see Equations 29-3 1). Without loss of generality, it has been assumed that the disc 
remains in its unrotated position (d = (/> = y = 0). 


When kQa»l, the function | S^(P{) appearing in Equation 1-D becomes very sharply peaked 
at the angle of reflection and the angle of transmission ( 0 p^(); and since W^^ (0,<f>) is a 

smoothly varying function under most conditions, the major contribution to for large k^a is 
expected to come from the vicinity of specular reflection and transmission. This leads to a pro- 
cedure to evaluate the integral in Equation 1-D which is asymptotically correct for large k^a. To 
formally carry out the procedure, the integral over the spherical surface surrounding the disk is 
projected onto two integrals in the plane of the disc by means of the transformation. 

= sin^cos 0 

5-D 

Oy = sin 0 sin 0 
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The transformed expression is: 






Ur At/ 


where 


6-D 


Af * jo^ + «y < I,cos5< o| 


Using the sharply peaked character of | S(Pp 1^, Equation 6D becomes 

CWW («„«,) (* !?(?;) P da,da^ t f | S(r;> |2 da^da^ 7-D 

Once again using the fact that kQa»I, the integrals in Equation 7D can be extended to the whole 
a„,a^ plane. The resulting integrals can then be evaluated exactly by employing Equation 22c. 

A y 

The result is 


where is the area of the disc. In arriving at Equation 6D the identity cos0j = I cos5^ ( = I cos^j 1 
has been used. 
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FIGURE CAPTIONS 


Figure 1. 
Figure 2. 
Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


Figure 8. 


Figure 9, 


Figure 10. 


Figure 1 1. 


Problem Geometry. 

Eulerian angles {6, <f>, y) which describe the disk orientation. 

Geometry as seen in the coordinate frame fixed on the disk. The origin is half-way 
between the parallel faces of the disk. 

Geometry used in calculating the examples, = 30*^ and 1 is in the y-z plane, n is 
also in the y-z plane and except where noted h > z as shoiVh here. 

Components of the scattering amplitude when the observer is in the plane of incidence 
(O'*), perpendicular to the plane of incidence (90®) and half-way between (45®). 

Magnitude of the scattering amplidue, in the plane of incidence as a function of 
ka. 

Magnitude of the scattering amplitude, f^y, in the plane of incidence as a function of 
ka. 

Magnitude of the scattering amplitude, as a function of disk orientation. 6 is the 
angle between n and i. 

Magnitude of the scattering amplitude, f^y, as a function of disk orientation. 6 is the 
angle between n and i. 

Magnitude of the scattering amplitude, fjjj^, in the plane of incidence as a function of 
dielectric constant (no loss): 7^ = e'. 

Magnitude of the scattering amplitude, fyy, in the plane of incidence as a function of 
dielectric constant (no loss): ej. = e'. 
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Figure 12. Magnitude of the scattering amplitude, f^|j, in the plane of incidence near resonance 
for dielectric with loss. 

Figure 12. Magnitude of the scattering amplitude, f^^, in the plane of incidence near resonance 
for dielectric with loss. 

Figure 14. Magnitude of the scattering amplitude, in the plane of incidence for purely lossy 
dielectric: 

Figure 15. Magnitude of the scattering amplitude, f^^, in the plane of incidence for purely lossy 
dielectric: 

Figure 16. Polarization and propagation vectors inside the slab. 
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